Comment on "Aging Effects in a Lennard- 
Jones Glass" 

In a recent Letter Kob and BarratS reported results of 
molecular dynamics simulations for the off-equilibrium 
dynamics in a binary Lennard- Jones (LJ) glass. The 
main conclusions of their work was 1) they find aging 
in this glassy systems and 2) that they find a simple ag- 
ing scenario close to a t/t w scaling, which is very rem- 
iniscent of comparable studies in spin glasses. In this 
comment we would like to emphasize that a different ag- 
ing scenario, known under the name activated dynamics 
scaling, is much more appropriate for the system under 
consideration than the one proposed by Kob and BarratE. 

For this reason we repeated the simulation by Kob and 
Barrat, using exactly the same potential (Lennard- Jones 
for a binary mixture), the same parameters (same diam- 
eters, mixture, density and temperatures) and the same 
quenching procedure (Tj=5, Tf=0A) however with much 
larger systems (32768 = 32 3 particles) and similar times 
(2 • 10 6 time steps, 1 time step corresponding to 0.01 
LJ-units). The aging properties of the system manifest 
themselves in the two-time autocorrelation function 
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where ri(t) is the position of particle i at time t and 
the absolute value of q corresponds to the first maximum 
in the structure function. We choose 100 randomly dis- 
tributed vectors and averaged C q over these vectors. The 
function ([!]) was evaluated after every 10 time steps and 
5™ measurements were averaged over to improve statis- 
tics. We convinced ourselves that different quenching 
procedures with identical initial and final temperatures, 
Tj and Tf, lead to the same scaling behavior. 

In [1] it has been suggested that C q (t + t Wl t w ) obeys 

C q (t + t w ,t w ) ~c(t/*r) (2) 

with a relaxation time t r cx i° . We checked this Ansatz 
for our data and display the result in the inset of Fig. 
1, surprisingly we find an exponent a ~ 1.1, very close 
to one (corresponding to simple t/t w scaling) but differ- 
ent from the one a = 0.88 reported in [1]. The data 
collapse in the asymptotic regime is not at all satisfying, 
the data for different waiting times conincide exactly only 
for C q = 0.45. For this reason we tried another aging sce- 
nario, proposed in the context of spin glasses by Fisher 
and Husecl, which we call the activated dynamics: 

C q {t + t w ,t w )~C{\n{{t + t w )/T)/Mt w /T)} (3) 
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FIG. 1. Activated dynamics scaling plot according to eq.(^) 
with r = 0.5. Note that we shifted the scaling variable by 1 
to the left to have a better resolution of the crossover region. 
The inset shows a scaling plot of our data according to the sce- 
nario proposed by Kob and Barrat [1], see eq.(^), the full line 
corresponds to C q = 0.45. The relaxation times are t r = 1450, 
10600, 70000, 600000 and 3000000 for i„=500, 2500, 12500, 
62500 and 312500, roughly a dependence t r oc . 

The origin of such an activated dynamics scaling in 
spin glass phenomenology!] is simply a logarithmically 
slow coarsening process £(t) ~ ln(t) a , where £(t) is a 
time dependent spatial correlation length and a some ex- 
ponent. This plus the observation that in coarsening dy- 
namics the two time correlation function C q (t + t w ,t w ) 
should depend on the ration of the two length scale 
£{tw)/£,(t + t w ) alone yields the aging behavior (||). 

Three things are worth being noted: 1) In the con- 
text, in which (||) was first suggested, nameljuthe 3d EA 
spin glass, this form does not seem to worm. 2) Only 
very recently a growing length scale has been observed 
in the very same model we are considering herea. 3) An 
even better data collapse can be obtained by plotting 
C q {t + t w , t w ) versus ln(i)/ ln(t r ), with a relaxation time 
t r individually chosen for each waiting time t w . Here it 
turns out that t r (t w ) grows faster than with a power law. 

To conclude we have shown that the aging behavior of 
a Lennard- Jones glass is more appropriately described by 
an activated dynamics scaling rather than simple aging, 
as claimed by Kob and Barrat in [1] . 
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where r is a fit-parameter and plays the role of an effec- 
tive microscopic time scale. Fig.l we show the scaling 
plot for such a scenario, which gives a much better data 
collapse in the asymptotic regime t > t w . 
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